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NOMENCLATURE 


Different  meanings  are  separated  by  semicolons. 

1  =  conductor  1,  forward  section,  also  used  as  a  subscript 

2  =  conductor  2,  rear  section,  also  used  as  a  subscript 

3  =  conductor  3,  extended  probe,  also  used  as  a  subscript 

a  =  radius  of  a  sphere  or  a  cylinder 


A,An  =  area,  area  of  conductor  n;  Ampere 


b  =  half  length  of  a  cylinder 

B  =  half  length  of  cylindrical  sheath 

c(a,b)  =  capacitance  of  a  conducting  right  circular  cylinder  a,b,  Eq.  (5) 

=  coefficient  of  capacity  or  induction,  clement  of  matrix  C,  Section  3 
=  capacity;  Coulomb 

=  3x3  capacitance  coefficient  matrix.  Section  3 
=  capacitance  of  the  forward  section  by  itself 


c.. 

iJ 


f 


=  capacitance  of  the  forward  section  with  the  rear  section  attached, 
insulated,  and  uncharged 


fr 


=  capacitance  of  the  forward  and  rear  sections  attached  and  electrically 
connected 


Vll 


kS 


m 


Sf 


•  %  * 
vMJ 
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Nomenclature  (cont'd) 
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D  =  P11P22  ~  p12  ’  determinant  P’  Section  3 

e  =  magnitude  of  electron  charge;  2.718...;  subscript  refers  to  electrons 

E  =  electric  field 

Eb  =  kinetic  energy  of  beam  particle 

e  =  electrical  conversion  factor,  1/(4tte  )  =  8.9876  m/nf 
6  o 

I  =  beam  current  leaving  forward  section 

I  =  net  (ion  and  electron)  current  from  plasma  to  conductor 

I  =  current  from  plasma  (excluding  beam  particles)  returning  to  conductor  n 

I  =  conventional  voltmeter  current  from  conductor  n  to  1 ,  Section  4 

nvm 

k  =  Boltzmann  constant  when  multiplied  by  T;  otherwise  1.777 
m  =  mass 

m  =  mass  of  an  electron 

e 

m.  =  mass  of  an  ion 

n  =  ambient  electron  or  ion  number  density 
o 

p„  =  coefficient  of  potential,  element  of  matrix  P,  Section  3 
p103  =  value  P13  would  have  if  conductor  2  were  missing 

p123  =  value  p13  would  have  if  conductor  1  and  2  were  connected 

P  =3x3  potential  coefficient  matrix.  Section  3 

q  =  electric  charge 

q.  =  charge  on  conductor  i,  element  of  Q;  instantaneous  charge  on  i,  Section 
1  4.10 

qg  =  charge  in  the  sheath 

Q  =3x1  column  charge  matrix 
r  =  spherical  radial  coordinate 

r  =  ratio  of  charge  on  conductor  1  to  that  on  2  when  both  are  at  the  same 
c  potential 

R  =  (usually  cylindrical)  radial  coordinate 


,  V  a-  -  ,’V 


Nomenclature  (cont’d) 


=  sheath  radius 


=  initial  ion  traverse  time  through  the  sheath.  Section  4.11 
=  potential;  Volt 

=  3x1  column  potential  matrix,  Section  3 
=  potential  of  conductor  i,  element  of  matrix  V 

=  3x1  column  matrix  of  conductor  potentials  due  to  sheath,  Section  3 
=  part  of  V.  that  is  due  to  the  sheath  charge,  element  of  V 
=  cylindrical  axial  coordinate 

=  resistance  of  voltmeter  between  conductors  i  and  1 

=  mean  number  of  secondary  electrons  per  primary  incident  particle, 
Section  4 

=  ratio  of  conductor  radius  to  Debye  length 
=  net  electrical  charge  density 

=  charging  time  constant  of  conductor  1  due  to  currents  and  capacitance, 
Section  4.8 

=  Vi  -  V1 
=  eV/kT 
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Theoretical  Analysis  of  Charging  Data 
From  Rocket  With  Charged  Beam  Emission 


1.  INTRODUCTION 


The  emission  of  an  ion  or  electron  beam  from  a  spacecraft  into  the  ambient 
plasma  changes  the  potential  on  the  beam-emitting  section  and,  by  induction,  on 
other  parts  of  the  spacecraft  relative  to  that  of  the  ambient  plasma.  The  values 
of  the  resulting  potentials  depend  on  the  parameters  of  the  beam,  rocket,  and 
plasma.  Since  this  dependence  was  not  well-known,  Cohen  et  al.  carried  out 
an  experiment  with  a  rocket  in  the  lower  F- region  of  the  ionosphere  at  night. 
Potential  differences  between  pairs  of  three  isolated  conductors  were  measured 
when  a  positive  ion  or  electron  beam  was  emitted  from  one  of  them.  These  volt¬ 
ages  seemed  remarkably  high,  even  for  relatively  low  positive  beam  currents. 
Questions  arose  as  to  why  some  of  the  measured  values  resulted.  This  report 

analyzes  some  of  the  data  to  try  to  answer  the  questions,  and  to  gain  further 

2  3 

understanding  of  the  results.  After  the  original  analysis,  Katz  and  Mandell 

(Received  for  publication,  October  198G) 

1.  Cohen,  H.  A.,  Sherman,  C.,  and  Mullen,  E.  G.  (1979)  Spacecraft  charging 

due  to  positive  ion  emission:  an  experimental  study,  G.R.L.  6:5-15. 

2.  Dubs,  C.  (12  Nov.  1982)  Potentials  and  Charges  on  Conducting  Rocket 

Sections,  AFGL-TR-82-0349 ,  ADA130143. 

3.  Katz,  I.,  and  Mandell,  M.  J.  (1982)  Differential  charging  of  high-voltage 

spacecraft:  the  equilibrium  potential  of  insulated  surfaces,  J.G.R., 
87:4533-4541. 


issued  a  paper  that  contains  some  analysis  of  this  experiment.  These  analyses 
are  compared.  This  report  is  a  revision  and  an  extension  of  the  original 
analysis. 

In  this  report,  "ion"  means  "positive  ion."  Section  2  outlines  the  experi¬ 
ment  and  lists  the  data  analyzed  here.  Values  of  coefficients  of  potential, 
capacity,  and  induction  are  determined  in  Section  3.  They  contribute  physical 
insight,  including  the  charge  on  each  conductor  and  in  the  sheath,  the  potential 
of  the  position  of  the  probe  if  the  probe  was  removed,  and  the  charging  time 
constant  due  to  capacitance  and  currents.  For  any  configuration ,  these  coeffi¬ 
cients  depend  only  on  the  geometry.  For  sufficiently  simple  geometry,  they  can 
be  calculated,  as  is  done  in  this  case,  without  direct  use  of  potential  theory. 

How  they  can  be  measured  before  or  after  flight  for  verification,  or  in  lieu  of 
calculation  for  any  geometry  and  number  of  conductors,  is  outlined.  In  Section 
4,  the  absolute  potentials  of  the  rocket  conductors,  charges  on  them,  sheath 
parameters,  probe  position  potential,  time  constants  for  steady  state,  etc.,  are 
determined  from  the  given  parameters.  Useful  empirical  formulae  are  also  given 
for  the  theories  of  Laframboise,  Langmuir-Blodgett ,  and  Lam.  Discussion  includes 
the  charge  and  potential  behavior  after  the  beam  is  turned  on,  comparison  of 
measured  currents  with  those  calculated  from  probe  theory,  secondary  electrons, 
sense  of  surface  charge  densities,  saturation  phenomena,  and  the  effect  of 
voltmeter  resistance.  The  conclusions  are  summarized  in  Section  5. 

2.  EXPERIMENT 

Cohen  et  al.1  flew  a  rocket  with  three  conductors  insulated  from  each  other: 

forward  section,  1,  rear  section,  2,  and  extended  thermal  emissive  probe,  3  (see 

Figure  1  and  Table  1}  to  258  km  altitude  at  night.  Conductors  1  and  2  were 

coaxial  cylinders  separated  by  an  insulating  ring.  Conductor  3  was  a  sphere 

and  a  short,  thin,  hot  wire  filament  (not  shown).  At  certain  times,  an  ion  or 

electron  beam  of  known  current,  1^,  and  particle  energy,  E^,  was  emitted  from 

the  front  of  the  forward  section,  causing  certain  steady-state  potentials  and 

charges  to  develop  on  the  three  conductors.  Voltmeters  measured  =  V2  -  V) 

and  ♦31  =  v3  -  vr  The  numbers  in  Table  2  are  the  part  of  the  results,  nine 

cases,  that  are  analyzed  here.  The  value  of  I,  designates  the  case.  The  plasma 
9-3  0 

density  was  near  10  m  and  the  temperature,  T,  of  both  ions  and  electrons  was 
near  550°K,  so  these  values  are  assumed  here.  The  potentials  were  measured 
after  steady  state  was  reached,  so  the  return  current  from  the  plasma,  minus 
secondary  electron  current,  plus  returning  beam  current  if  above  saturation 


A 
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Table  1.  Rocket  Dimensions  in  Meter 
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a3 

bl 

b2 

b12 
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.19 

.19 

.019 

1.27 

.205 

1.475 

.26 

1.71 

Table  2.  Measured  Beam  Current  and  Energy,  Rocket  Altitude,  and  Conductor 
Potential  Differences 
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Elect  i  on 


ohen  Table 


Fig.  Table  Fig 
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(Section  4),  equals  Ij  ,  and  electrostatic  theory  is  used,  except  for  charging  and 
ion  transit  times.  An  ion  beam  was  emitted  with  the  electron  beam,  but  its  cur¬ 
rent  was  negligibly  small;  thus,  it  is  ignored  in  the  "electron"  column  of  Table  2. 
At  least  for  ion-beam  emission,  t lie  data  were  found  to  be  independent  of  pitch 
angle  and  neutral  number  density.  For  this  reason,  as  well  as  for  simplicity,  the 
magnetic  field  and  neutral  density  arc  ignored  in  the  analysis.  The  effect  of  ram 


is  neglected  since  the  rocket  orientation  is  unknown,  but  it  should  be  important 
for  ion  current.  Questions  arose:  What  were  the  absolute  values  of  the  poten¬ 
tials?  Why  was  4^  <  The  purpose  of  this  report  is  to  make  an  approxi¬ 

mate  analysis  of  the  measurements  in  Table  2,  and  to  gain  an  understanding  of 
some  of  the  basic  electrical  phenomena,  as  well  as  to  try  to  answer  these 
questions . 


3.  COEFFICIENTS  OF  POTENTIAL,  CAPACITY,  AND  INDUCTION 


The  purpose  of  this  section  is  to  determine  values  of  the  coefficients  of 

potential,  capacity,  and  induction  for  the  three  conductors  of  the  rocket  system. 

These  coefficients  relate  the  charges  to  the  potentials  on  these  conductors.  Since 

they  were  not  measured,  they  are  calculated.  Fortunately,  the  geometry  is  simple 

enough  so  that  it  is  not  necessary  to  start  from  potential  theory.  The  coefficients 

4 

are  obtained  by  using  Smythe's  formula  for  the  capacitance  of  a  conducting  right 
circular  cylinder,  an  empirical  formula  found  for  the  potential  at  a  point  near  such 
a  charged  cylinder,  carrying  out  three  gedanken  experiments,  and  interpolating. 

Let  p^,  c..,  and  c^  for  i  *  j  be  the  coefficients  of  potential,  capacity,  and 
induction,  respectively,  the  elements  of  matrices  P  and  C.  From  electrostatics, 
e.g.,  a  slight  generalization  of  Page,5 


V  =  PQ  +  Vs. 

(1) 

Q  =  C(V-VS). 

(2) 

C  =  p”1, 

(3) 

P  and  C  are  symmetric,  the  p.^  and  c./s  are  positive,  and  the  c..'s  for  i  *  j  are 
negative.  Vg  accounts  for  the  effect  of  the  plasma. 

Physical  reason  leads  to  the  inequality  of  rocket  capacitances 


<  C 


fr ' 


(4) 


4 

Smythc  gave  the  following  empirical  formula  (here  slightly  modified  to  make  it 
correct  within  .30%  instead  of  .41%)  for  the  capacitance  of  a  conducting  right 
circular  cylindrical  conductor  of  radius  a  and  half-length  b  for  a/8  ■_  b  *  8a: 


4.  Smythe,  W.  R.  (1962)  Charged  right  circular  cylinder,  J.  Appl.  Phys. 
33:2966-2967. 


5.  Page,  L.  (1935)  Introduction  to  Theoretical  Physics,  Second  Edition,  §  120. 


c(a,b)  =  [.0707  +  ,0615(b/a)  ]a  nf,  (5) 

a  and  b  in  meter.  In  particular,  from  this  and  Table  1, 

Cf  =  .0629  nf  (6) 

Cfr  =  .0689  nf.  (7) 

Three  gedanken  experiments,  designated  with  the  subscripts  a,  b,  and  c, 
are  now  carried  out  in  a  vacuum.  Eq.  (1)  with  Vg  =  0  is  used  for  each.  In  the 
first  experiment,  charge  is  put  on  1  while  2  and  3  are  left  uncharged.  Then 

Vla  =  Pll^la  (8) 

V2a  =  Pl2^1a  (9) 

V3a  =  Pl3^1a'  (10) 

When  uncharged,  3  is  small  and  far  away  enough  from  1  and  2  to  affect  their 
potentials  negligibly,  so 


1l.  -  Cr'via' 


■>11  =  1/cf'-  02 

Physical  consideration  shows  that  |V2a]  <  S0  ^s'  ^  an( d  (9)  show  that 

P12  <  Pir  (13 

In  the  second  experiment,  charge  qlfa  is  put  on  1  and  q2b  on  2  such  that  Vlb  - 
V2b  with  no  charge  on  3.  Then 

Vlb  =  Pll^lb  +  P12q2b  (14 

Vlb  =  PlZ^lb  +  P22q2b  (15 

V3b  =  Pl3^1b  +  P23q2b-  (16 

Again,  3  is  small  and  far  away  enough  so  that 

Vlb  =  'lib  *  <WCfr"1  (17 

accurately.  This  and  qJb  =  rcq2b  are  substituted  into  EQS>  d4>  and  d5). 

The  former  may  then  be  written 

Pl2  =  (rc  *  1  )Cfr_1  '  rcpir  (18 


8 


| 

1 

I 


Substituting  this  into  the  latter  leads  to 

p22  =  rc2pll  "  (rc2  "  1)Cfr_1*  (19) 

4 

How  is  rc  obtained?  Smythe  expressed  the  surface-charge  density  of  a  charged 
conducting  right  circular  cylinder  as  fractional  power  series,  in  z  for  the  side, 
and  in  R  for  the  ends  (inner  cylinder.  Figure  2).  He  obtained  the  coefficients 
for  several  values  of  b/a.  The  coefficients  for  this  ratio  equals  8  are  used  since 
t*12^al  =  7.76.  and  q2  are  obtained  by  integrating  these  surface  densities 
over  1  and  2  separately.  The  ratio  gives 

rc  =  3.770  (20) 

Details  of  this  calculation  are  in  Reference  2,  Section  Al.  Thus,  p^,  p^,  and 
P20  may  be  obtained  given  C^'. 

From  Eq.  (10), 

Ptf'*  -  <lla/V3a-  <2l> 

The  value  p13  1  would  have  if  2  were  missing  is  shown  in  Appendix  A  and  is 
given  by  Eq.  (A5), 

Pl03_1  -  •2605  nf'  <22> 

Since  3  is  small  and  far  enough  away,  if  1  and  2  were  electrically  connected, 
p,,  1  would  be 


_  -1  _  qlb  +  q2b  , 

p‘22 


Eq.  (A6),  Appendix  A.  Physical  reason  shows  that 


-1  n  -1  -1 
P103  P13  P 1 23  ’ 


similar  to  Eq.  (4).  from  Eq.  (23)  and  q^  =  r^q,^  are  substituted  into 

Eq.  (16),  yielding 

p23  =  (rc  +  1)p123  '  rcp13‘  ( 


m 

m 

sty 

:$ 


i 
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In  the  third  experiment,  charge  q3c  is  put  on  3  with  no  charge  on  1  or  2. 


Then 


V3c  p33q3c‘ 


A  maximum  effect  calculation  shows  that  the  presence  of  1  and  2  changes  V3c  by 
less  than  .1%,  so  this  effect  is  neglected.  Therefore,  q3c  =  C3V3c,  where  the 
capacitance  of  3  is  =  a3/g  +  C  ,  where  is  the  capacitance  of  the  hot  wire. 


For  the  wire,  2a  =  .00381  cm,  2b  =  3.175  cm,  so  C  =  .0002620  nf,  obtained  for 


q/V  from  Eq.  (Al)  with  R  =  a  and  z  =  0,  Appendix  A.  (It  could  also  be  obtained 
from  the  formula  for  the  capacitance  of  a  wire,  the  diameter  of  which,  2a,  is  small 
compared  to  its  length,  b/[g  ln(b/a)]  =  .0002625  nf.  Smythe's  formula,  Eq.  (5), 
cannot  be  used  since  it  does  not  reduce  to  the  latter  for  b  >>  a).  This  is  11%  of 


C^.  So,  from  the  value  of  a^  in  Table  1, 


P33  =  420  (nf) 


Thus,  all  p.-'s  may  be  obtained  given  C^.'  and  p^3 

With  Eqs.  (4)  and  (24)  in  mind,  the  value  b.  +  b„/2  of  effective  half-length 

1  i  _i 

of  conductor  1  due  to  2  being  present,  and  the  average  of  plfl3  and  P123  f°r 
p,,  1  are  chosen.  Then 


Cf'  =  c ( a j ,  bj  +  b2/2)  =  .  0659,.  nf. 


P13  =  .261  nf. 


These  result  from  using  the  values  in  Table  1,  Eqs.  (5),  (22),  and  (23).  The 
value  in  Eq.  (28)  accurately  equals  the  average  of  the  values  in  Eqs.  (6)  and 
(7).  From  the  above, 


12.1  3.83 


23.7  3.79 


3.79  420 


These  values  should  be  correct  within  5%,  assuming  negligible  errors  in  the 
numbers  in  Table  1.  From  Eqs.  (3)  and  (30), 


111. 

-56.2 


-56.2  -.503 

70.7  -.126 


Pf. 


(31) 


Ignoring  3,  a  good  approximation, 

C11  =  p22/D’ 
c12  =  C21  =  ~P12/D’ 

C22  "  Pll/D’ 

D  =  P11P22  -  P122  =  <rc  *  "V*  -  Cfr_1] ' 

Choosing  q„  -  -  q^,  the  capacity  between  1  and  2  is 


(32) 

(33) 

(34) 

(35) 


ql^Vl  V2^  1  ^ p  1 1  ”  2p12  p22^ 

=  (C11C22  '  C122)/(C11  +  2C12  +  C22)  =  1/(CfrD)‘  (36) 

So,  as  the  gap  between  1  and  2  approaches  0,  tins  capacity  and  | c. ^  1  become 
very  large  since  Cy  approaches  C^.  Including  3,  as  Cj.'  approaches  C^,  D 
approaches  0  and  pJ3  approaches  p^. 

Since  the  value  of  C^.'  can  be  critical,  the  question  naturally  arises,  "How 
reasonable  is  Eq.  (28)?"  Rocket  sections  1  and  2  were  separated  by  a  2-inch 
thick  cylinder  of  fiberglass  that  was  essentially  hollow  except  for  a  one-half-inch 
thick  fiberglass  disc.  To  check,  the  capacitance  between  the  two  ends  of  an 

18.7  cm  length  of  the  rocket  that  included  the  gap  was  measured  and  found  to  be 

44.7  pf.  A  correction  due  to  the  missing  lengths  of  1  and  2  is  calculated  to  be 
10.2  pf.  Adding  this  and  equating  to  Eq.  (36),  D  is  obtained,  which  from  Eq. 
(35)  leads  to  Cj.'  =  .0653  nf,  insensitive  to  the  value  of  capacity  between  1  and  2, 
close  to  the  value  in  Eq.  (28). 

These  calculations  of  P  can  be  checked  or  replaced  by  six  measurements  of 
capacity  by  hanging  the  rocket  over  10  times  its  largest  dimension  (for  <  5?, 
error)  from  comparable  and  larger  size  conductors  and  ground,  and  running  a 
small  ground  wire,  0,  to  the  vicinity  of  the  rocket.  Measurement  of  the  capa¬ 
cities  between  1  and  0,  2  and  0,  and  3  and  0  give  the  p..'s.  Measurement  of  the 
capacities  between  1  and  2,  2  and  3,  and  3  and  1  give  the  p.?s  for  i  *  j. 

How  this  may  be  generalized  to  any  number  of  conductors  is  evident . 
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4.  NUMERICAL  RESULTS 


This  section  contains  calculations  of  numerical  values  of  at  least  ten  types  of 
physical  quantities  and  a  discussion  of  them.  Unlike  Reference  2,  probe  theory 


is  used.  Where  applicable,  calculated  results  are  compared  with  experiments,  and 

3 

deductions  are  compared  to  those  of  Katz  and  Mandell. 


4.1  .Absolute  Potentials 


First,  the  potentials,  V.,  of  the  three  conductors  relative  to  that  of  the 
ambient  plasma  are  calculated.  This  is  done  by  determining  the  values  of  V2> 


then  using  the  measured  values  of  4.^  from  Table  2.  In  Reference  2,  it  was 


assumed  that  V2  =  0.  V2  must  be  sufficiently  negative  to  attract  a  return 
current  equal  to  the  voltmeter  current,  I„  =  $21^Z12’  P*us  net  current  due  to 


secondary  electrons.  The  latter  current  is  neglected,  but  it  could  be  important. 


z12  =  io  a. 


The  formulas  used  here  that  relate  the  return  current,  1,  due  to  ions,  I., 

i 


and  to  electrons,  Ig,  from  the  surrounding  plasma  to  a  conductor  at  a  negative 


potential  V,  are  the  following: 


x(L) 


- 


=  e 


(37) 

(38) 


(39) 


I  =  I.  +  I 
l  e 


I.  =  I.  i. 

i  io  i 


I  =  I  i 

e  eoe 


I.  =  n  ev..A 
io  o  ti 


(40) 

(41) 

(42) 

(43) 


eo 


Vm. 

—  I. 
m  io 
e 


(44) 


kT 


(45) 
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S3 


Eq.  (39),  (40),  (41),  (42),  or  (44)  is  redundant.  Tg  =  =  T  is  assumed  since 

ion  and  electron  temperatures  are  nearly  the  same  in  this  experiment. 

Since  Table  1  shows  a2  and  b2  to  be  nearly  equal,  spherical  probe  theory  is 
used  for  Eq.  (37)  to  determine  V2.  The  effect  of  V1  on  I2>  the  current  to  2,  is 
neglected.  Probably  being  negative  only  cuts  down  the  electron  collecting 
area,  affecting  only  the  1  uA  case.  For  1^  =  1  iiA  or  -10  mA,  [V2|  is  small 

enough  for  the  results  of  Laframboise®  to  apply.  For  £  <  5,  the  following  fit  to 

6 

Table  5c  of  Laframboise  is  used: 


1  -  1<1 

=  kjdj  -  1)  +  — -  ln[  1  +  k2(i.  -  1)] 

2 

(47) 

,  .  6.836  -  .685£2 

1  - n — 

(5  -  6) 

(48) 

=  .00516(5  -  6)1'7  • 

(49) 

These  have  an  error  of  less  than  1%  for  £  <  1,  2%  for  6  <  3,  and  4%  for  £  <  5. 
Interpolated  correction  factors  are  used  so  that  the  error  is  cut  to  -  .2%.  A 
more  accurate  interpolation  and  extrapolation  for  £  <  .5  is  given  in  Section  4.9. 

1  =  I  ,  and  the  conductor  area,  A,  and  radius,  r„,  are  chosen  to  satisfy  A  = 

2  i  2  i  ,  i  4 

4rrr2  =  4-a2b2  +  -a2.  £  =  4.28.  The  results  are  |  V g !  in  the  lb  =  1  and  -10  uA 

columns  of  Table  3.  For  each  of  the  other  values  of  beam  current,  |V’2j  is  so 

large  that  i  is  negligible,  so  is  set  equal  to  zero  in  Eqs.  (38)  to  (42),  and  x  is 

e  7 

out  of  the  range  of  Laframboise's  results.  For  them,  Lam's  "highly  negative" 

g 

spherical  probe  potentiiil  theory  is  used  for  Eq.  (37).  (The  Langmuir-Blodgett 

1/3 

formula  for  V  is  (4/-)  times  the  Lam  formula  with  the  argument  of  F  multiplied 

2 

by  1.26  (larger  unless  £  j  >>  16),  so  it  gives  larger  values  for  'V2I.  It  is  not 
used  since  it  disagrees  much  more  with  values  extrapolated  from  Laframboise's 
results,  and  since  it  leads  to  currents  further  from  experimental  values.) 


6.  Laframboise,  J.  G.  (1966)  Theory  of  spherical  and  cylindrical  Langmuir 

probes  in  a  collisionless,  Maxwellian  plasma  at  rest,  UT1AS  Report  No.  100. 

7.  Lam,  S.  H.  (1965)  Unified  theory  for  the  Langmuir  probe  in  a  collisionless 

plasma,  Phys.  Fluids  8:73-87. 


8.  Langmuir,  I.  and  Blodgett,  K.  B.  (1924)  Currents  limited  by  space  charge 
between  concentric  spheres,  Phys.  Rev.  24:49-59. 


Given  F,  e  from  Eq.  (54)  is  correct  within  2/3%  for  0  <  z  <  1.  The  resulting 

values  of  V„  for  I.  =  8,  9,  and  12  uA  are  in  Table  3.  How  the  values  are  ob- 
l  b 

tained  for  =  10.9  and  374  uA  is  explained  later.  The  values  of  and  in 
Table  3  follow  from  Table  2. 

The  values  of  V2  are  seen  to  be  relatively  close  to  zero.  This  is  to  be 
expected  since  Z21  is  large.  After  an  ion  (electron)  beam  is  turned  on,  the 
positive  (negative)  charge  removal  causes  an  increasing  negative  (positive) 
charge  on  1  that  induces  a  negative  (positive)  potential  on  2,  Thus,  both  sec¬ 
tions  attract  return  ions  (electrons)  from  the  plasma.  In  steady  state,  1  must 
be  sufficiently  negative  (positive)  to  attract  a  return  current  equal  to  the  beam 
current  minus  the  secondary  electron  current,  here  neglected,  minus  I9.  Conduc¬ 
tor  2,  however,  must  be  only  negative  enough  to  attract  a  return  current  to  can¬ 
cel  voltmeter  and  secondary  electron  current  from  1  to  2.  The  secondary  electron 
current,  here  neglected,  effectively  decreases  Z21  making  V2  more  negative.  So, 
qj  and  q2  must  induce  nearly  equal  but  opposite  potentials  on  2.  The  sign  of  q2 
is  the  same  as  that  of  the  beam  particles,  opposite  to  that  of  q^,  and  | q 2 j  <  Iq^. 
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4.2  Return  Current  to  Forward  Section 

Lam's'  "highly  negative"  cylindrical  probe  potential  theory  is  used  to 
calculate  the  return  current  from  the  plasma  to  1,  Ij,  =  I^am’  ^rom  eac^  value  of 
in  Table  3. 


*1 


(56) 


V 


4  +  ir.943 


‘10 


1 


(57) 


G(x)  = 


2/3 


(58) 


2  9 

where  (-6)  as  a  function  of  x  =  rQ/r  is  given  by  Langmuir  and  Blodgett.  An 

empirical  formula  found  for  G  correct  to  1%  is 


(1.5y)4/3  +  .869y3 
G(x)  =  - r 

1  +  . 271y  +  .417y  (59) 


m 

i 


m 

■rfyn 


C:>i 


a 


V  V 


■f 


Y  =  In  x.  (60) 

9 

(This  y  is  the  negative  of  Langmuir- Blodgett's  y)-  Given  G,  x  -  1  from  Eq. 

(59)  is  correct  within  3/4%  for  1  <  x  <  2.  The  values  of  I,,  obtained  are  in 

lLam 

Table  3. 

For  the  1^  =  10.9  and  374  yA  cases,  was  unknown,  so,  by  starting  with 
an  initial  intelligent  guess  and  iterating,  1^,  I2,  V^,  and  V2  are  determined  so  as 
to  satisfy  x^d^)  from  Lam's  cylindrical  probe  theory,  X2(^  ^rom  Lam's  spherical 

probe  theory,  and  $2^  =  ^^12'  ^or  *b  =  ^b  ~  ^1  +  *2  *s  8480  used*  For 

Ib  =  374  y  A ,  the  trial  assumption  is  made  that  the  case  is  saturated  (e|V^|  = 

E^),  so  Vj  =  -2000  V.  The  value  of  follows  from  and  it>^ ^  from  Table  2. 

The  values  obtained  for  these  five  quantities,  for  each  of  these  two  cases,  are 
shown  in  Table  3.  For  1^  =  10.9  yA,  is  out  of  line  with  other  values  of 
because  #3^  is  out  of  line. 

9.  Langmuir,  I.  and  Blodgett,  K.  B.  (1923)  Currents  limited  by  space  charge 
between  coaxial  cylinders,  Phys.  Rev.  22:347-356. 
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For  beam  currents  from  1  through  12  u A ,  1 


is  seen  to  agree  roughly 


with  Ib  -  1 2 .  so  Lam  theory  accounts  for  the  high  negative  potentials  of  1.  The 
discrepancy  for  low  values  of  1^  could  be  due  to  simplifying  assumptions.  Ions 
incident  on  the  forward  section  have  an  energy  -eV^.  At  100  km  altitude  and  33n 
latitude,  returing  NO+  ions  moving  perpendicularly  to  the  magnetic  field  have  a 
gyro-radius  of  450  (150)  m  at  500  (50)  eV  energy.  For  other  angles,  these  radii 
are  decreased  by  the  factor  sin  u,  where  u  is  the  pitch  angle.  Since  the  lower 
energy  ions  are  more  inhibited  from  returning  by  the  magnetic  field,  perhaps  this 
discrepancy  is  due  to  the  geomagnetic  field.  Alternatively,  the  measurement  of 
may  have  been  low  at  the  small  values,  or  a  combination  of  both  of  these  factors 
could  have  caused  the  discrepancy. 

4 

For  beam  measurements  of  374  anti  -10  aA,  's  socn  *°  ^e  muc^  l°ss 

than  I  -  I2>  so  these  cases  were  highly  saturated;  most  of  the  returning  current 
is  beam  current  reflected  back  near  the  sheath  boundary.  For  the  374  ^A  case, 

3 

this  is  consistent  with  the  interpretation  of  Katz  and  Mandell,  Section  4.1. 

I. 3  Sr<:on(lar\  Electrons  and  Ionization 

The  effect  of  secondary  electrons  is  now  considered  briefly.  The  rocket 

surface  was  iridite-treated  aluminum.  For  the  ion-beam  cases,  values  of  ;,  the 

coefficient  of  secondary  emission  of  electrons,  for  NO+  ions  on  this  surface  for 

energies  up  to  2  keV  are  needed.  In  lieu  of  this,  the  value  assumed  is  .18,  the 

yield  from  1  keV  G+  ions  on  molybdenum.*®  All  of  the  secondary  electrons  from  1 

escape,  so,  Ijj  should  have  been  compared  to  (1^  -  I.,)/(l  +  1 ) .  This  slightly 

increases  the  discrepancy  at  small  values  of  I  Assuming  *  r  .18.  the  secondary 

electron  density  at  the  surface  of  1  is  -  17,  of  the  ion  density,  less  away  from  the 

surface,  so  conductor  and  sheath  charges  are  affected  negligibly. 

The  electron-beam  case  is  much  different.  Using  the  formula  of  Sternglass 

and  the  values  for  aluminum  :  •  -  .97  and  U  300  V,  •  -  .72  is  obtained 

max  max 

for  U  -  90  V.  Since  the  secondaries  have  but  a  few  cV  energy,  none  escape,  so 
the  return  current  is  unchanged.  If  all  of  the  secondaries  left  the  surface  with 
3  eV  of  energy,  the  electron  density  at  the  surface  would  be  nine  times  as  great 
as  with  6  =  0,  and  the  secondaries  would  extend  -  1/4  m  from  t lie  surface. 

10.  Stannard,  P.  R.,  Katz,  I.,  Mandell,  M.  J.,  Cassidy,  J.  J.,  Parks,  D.  E., 
Rotenberg,  M.,  and  Steen,  P.  G.  (1980)  Analysis  of  the  Charging  of  t lie 
SCATHA  (P78-2)  Satellite,  Report  NASA  CR-165348,  SSS-R-81-4798, 
pp.  22-24,  29. 

II.  Sternglass,  E.  J.  (1954)  Backscattering  of  kilovolt  electrons  from  solids, 

Phys.  Rev.  95:345. 
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Nevertheless,  the  magnitude  of  the  total  charge  of  the  secondaries  is  small  com¬ 
pared  to  |q  |,  so  and  qg  are  little  affected.  Some  incident  electrons  are 
backscattered.  The  number  is  less  than  the  number  of  secondaries,  and  they  are 
finally  captured. 

Ionization  by  beam  and  return  ions  is  less  than  a  few  percent  of  the  peak 
value  due  to  small  ion  speeds,  so  is  negligible.  For  the  electron-beam  case, 
however,  the  outward  and  returning  beam  and  the  return  electrons  from  the 
plasma  produce  ionization  in  the  sheath,  especially  close  to  the  forward  section. 
The  minimum  ionization  mean  free  path  (at  the  peak  of  the  ionization  cross  section 


vs  energy  curve)  is  400  m  at  138  km,  and  1300  m  at  160  km  altitude;  so,  it  is 
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expected  to  have  a  moderate  effect  at  the  most.  Since  ionization  constitutes  an 


additional  return  current,  the  sheath  radius  must  be  diminished  to  maintain  the 
same  total  return  current.  The  effect  is  not  treated  quantitatively  here.  The 
magnetic  field  causes  the  electrons  to  gyrate,  and  thus  inhibits  their  flow  across 
the  field.  This  effect  is  not  treated  quantitatively  here  either. 


4.4  Discussion 


The  values  of  V 2  in  Table  3  for  ion-beam  emission  are  relatively  small  and 
negative  in  qualitative  agreement  with  the  value  in  Case  I,  Table  2,  of  Katz  and 

3 

Mandell.  Except  for  1^  =  1  uA,  however,  they  are  appreciably  more  negative 


than  indicated  by  their  value.  This  is  due  mainly  to  taking  the  voltmeter  cur¬ 
rent  into  account  here.  The  front  part  of  2  had  a  positive  surface  charge  den¬ 
sity.  It  received  most  or  all  of  the  secondary  electrons  that  went  from  1  to  2. 

If  10%  (probably  much  too  high)  of  the  secondaries  from  1  went  to  2,  Ir  would  be 
increased  19  to  33%  by  increasing  A  from  0  to  .18.  The  sign  of  the  charge  on  the 
rear  part  of  2  may  be  determined  by  comparing  the  charge  on  2  that  is  tied  up  in 
the  capacitance  between  the  two  sections,  with  the  total  charge  on  2.  The  former 
is  the  product  of  and  this  capacitance.  (A  value  of  the  latter  consistent  with 
Cy,  Eqs.  (35)  and  (36),  p./s,  and  q.'s  should  be  used.)  The  resultant  product 
is  1.3  to  1.5  times  q„.  So,  the  charge  on  the  rear  part  of  2  is  concluded  to  be 

1  3 
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negative,  contrary  to  the  statement  of  Katz  and  Tlandell,  Section  3.1.  Trajec¬ 
tories  of  plasma  ions  existed  to  all  parts  of  2  except  possibly  to  the  (unknown) 
wake  side.  The  (— . 3Ig )  plasma  electron  current  for  the  1  yA  case  must  have 
reached  mainly  the  negative  surface  charge  density  part  of  2  because  of  the 
negative  potential  barrier  near  the  other  parts.  For  the  -10  mA  case,  the 


12.  Leadon,  R.  E. ,  Woods,  A.  J.,  Wenaas,  E.  P.,  and  Klein,  H.  H.  (1981) 
Analytical  Investigation  of  Emitting  Probes  in  an  Ionized  Plasma, 
AFGL-TR-81-01387  ADA104166,  §  3.4.2. 
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forward  and  rear  section  potentials  were  of  opposite  signs,  so  there  was  a  zero- 
potential  surface  surrounding  2  that  closed  on  the  separating  insulator.  As  with 
ion-beam  cases,  but  with  opposite  signs,  the  front  part  of  2  had  a  negative  sur¬ 
face  charge  density  and  the  rear  part  a  positive  surface  charge  density.  Second¬ 
ary  electrons  from  1  could  not  escape,  and  those  from  2  were  negligible  since  V2 
was  so  close  to  zero.  The  ionization  of  neutral  air  molecules  by  beam  and  return¬ 
ing  electrons,  here  neglected,  similar  to  increasing  nQ,  is  expected  to  increase 

and  |q  |  significantly,  but  not  to  affect  the  conductor  potentials  appreciably, 
s 


4.5  Sheath  Radii 


Sheath  radii  are  calculated  by  assuming  that  the  plasma  return  current 
equals  the  thermal  current  density  of  the  attracted  specie  times  the  sheath  area. 
When  the  sheath  is  small,  it  is  cylindrical  (Figure  2),  so 


noevt(2*Rsc2B  +  »Rsc>  =  II’ 


where 


The  latter  is  the  simplest  expression  found  for  the  half-length  of  a  cylinder  that 


increases  from  b  to  R  with  dB/dR  non-negative  and  increasing  monotonically 

SC  sc 


from  <  1  to  1  as  R  increases  from  a  to  R  >>  b.  When  the  sheath  is  large 

cp  —  cr*  ° 


compared  to  the  rocket,  the  sheath  is  spherical  (Figure  3),  so 


n  ev.4irR  =  I,, 
o  t  ss  1 


Except  for  1^  =  10.9,  374,  and  -10  yA,  these  sheath  radii  are  calculated  both  for 


1=1.-  I0:R  &  R 
1  b  2  scm  ssm 


,  and  for  1.  =  I1T  :  R  .  &  R  , 

1  lLam  sc  Lam  ssLum 


The  results 


are  shown  in  Table  4.  Column  1  ,  and  rows  R  and  R  ,  are  for  1,  =  I.  -  In. 

m  scm  ssm  1  b  2 


For  the  ion-beam  cases,  these  would  be  8  to  111  larger  if  6  were  assumed  to  be 
.18  instead  of  0.  For  the  electron-beam  case,  they  are  unaffected  by  6.  With 
the  possible  exception  of  the  1^  =  1  yA  case,  the  probe  (at  R  =  1.71  m)  is  seen 
to  be  inside  the  sheath.  This  agrees  with  the  possible  interpretation  in 
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Rtf 


io 


Rtf 


& 


Figure  3.  Cylinder  With  Spherical  Sheath 
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Table  4.  Sheath  Radii  and  Sheath  Potential  Coefficients 


!b 

1 

m 

1 

8a 

8b 

8c 

9 

10.9 

12 

374 

1 

O 

4* 

ua 

*1 

.95 

2.58 

11.7 

9.6 

10.7 

9.5 

10.4 

11.3 

25.8 

-853 

pa 

* 

scm 

1.3 

4.1 

4.2 

4.1 

4.4 

5.2 

m 

RscLam 

2.3 

5.2 

4.7 

5.0 

4.7 

4.9 

5.1 

7.9 

2.8 

m 

*ssm 

1.7 

4.9 

4.9 

4.9 

5.2 

6.0 

m 

RssLam 

2.9 

6.1 

5.5 

5.9 

5.5 

5.8 

6.0 

9.1 

3.4 

m 

plsc 

6.8 

4.3 

2.1 

2.3 

2.2 

2.3 

2.2 

2.1 

1.4 

3.7 

V/nC 

p3sc 

4.0 

3.6 

2.0 

2.2 

2.1 

2.2 

2.1 

2.0 

1.4 

3.3 

V/nC 

Plss 

7.6 

4.7 

2.2 

2.4 

2.3 

2.4 

2.3 

2.2 

1.5 

3.9 

V/nC 

P3ss 

5.3 

4.1 

2.1 

2.4 

2.2 

2.4 

2.3 

2.2 

1.5 

3.6 

V/nC 

pls 

7.2 

4.5 

2.2 

2.4 

2.3 

2.4 

2.3 

2.2 

1.5 

3.8 

V/nC 

p3s 

4.7 

4.0 

2.1 

2.3 

2.2 

2.3 

2.2 

2.2 

1.5 

3.5 

V/nC 

Except  for  the  1  column  and  the  I.  ,  R  ,  and  R  rows,  the  numbers  are  for 
•  _  m  b  scm  ssm 

*1  "  ‘lLarrT 


Section  4.1  of  Katz  and  Mandell.  For  the  1.  =  1  pA  case,  other  more  realistic 

b 

sheath  models  are  tried.  The  distance  of  the  sheath  edge  from  the  rocket  axis 
along  the  line  to  the  probe  is  calculated  assuming  1^  =  2.58  pA  and  the  sheath 
area  equals  Ij/nQevt.  For  the  forward  section  concentric  with  a  prolate  spheroid 


sheath  having  semi-minor  and  semi-major  axes  a  and  bQ  with 


b  =  a 
e  e 


2.6  m  is  obtained.  The  major  and  rocket  axes  are  coincident;  aQ  =  2.7,  bg  = 
3.2  m.  With  the  same  sheath,  but  with  the  forward  section  displaced  along  its 
axis  so  that  its  rear  edge  is  a  distance  of  .41  and  .025  m  from  the  sheath  edge 
parallel  to  the  axis,  1.6  and  1.1  m,  respectively,  are  obtained.  The  last  model 
tried  is 


r  =  3.03  +  1.74cos  6  (65) 

in  a  plane  containing  the  axis  rotated  about  the  axis,  the  origin  at  the  center  of 
the  forward  section,  and  9  measured  from  the  forward  sense  of  the  axis.  For  y  = 
.025  m,  2.0  m  is  obtained.  Thus,  the  probe  is  concluded  to  be  at  or  near  the 
sheath  edge.  If  1^  =  2.58  pA,  the  results  of  the  different  models  differ  too  much 
to  conclude  whether  or  not  the  probe  was  inside  or  outside.  If  1^  =  .95  pA,  the 
probe  was  almost  certainly  outside  of  the  sheath. 

4.6  Sheath  Potential  Coefficients 

Values  of  potential  per  unit  sheath  charge  due  only  to  sheath  charge, 
assuming  a  cylindrical  sheath,  p.  ,  and  assuming  a  spherical  sheath,  p.  ,  are 

JSC  Job 

calculated  at  the  surface  of  the  forward  section,  j  =  1,  and  the  probe,  j  =  3,  as 

indicated  in  Appendix  B.  They  are  shown  in  Table  4  and  are  a  function  only  of 

the  geometry,  including  that  of  the  sheath.  The  mean  of  p.  and  p.  is 

JSC  JSS 

calculated,  weighted  according  to  sheath  radius: 

Pjs  =  (1  '  F)Pjsc  +  Fpjss  (66) 

F  _  Rsc-  al 

Rsc  -  bL  -  28l  .  (67) 

The  constants  in  Eq.  (67)  are  chosen  so  that  F  =  .5  when  R  =  b^.  Values  of 
p.  ,  j  =  1,3,  are  shown  in  Table  4.  For  I.  =  8,  9,  and  12  pA,  values  of  p.  for 

JS  D  JS 

lj  =  lb  -  lg  differ  from  those  for  1^  =  IjLam  by  half  t,ie  Percentage  that  the 
latter  current  differs  from  the  former  (Table  3). 


WWW 


4.7  Conductor  and  Sheath  Charges 


The  charges  on  the  conductors  and  in  the  sheath  are  now  obtained  for  each 
case.  In  steady  state,  the  total  rocket  and  sheath  charge  vanish;  therefore,  the 
sheath  charge  is 


%  =  "  (<*1  +  <*2  +  V- 

The  potential  of  the  conductors  due  only  to  sheath  charge  is 

VjS  =  pjsV 

where  the  approximation  is  made  that 
p2s  =  pls- 

From  these  equations  and  Eq.  (1), 

V  =  (P  -  Pg)Q, 

where 


(68) 


(69) 


(70) 


(71) 


=  plg  for  i  =  1  or  2 

=  p3g  for  i  =  3.  (72) 

So, 

Q  =  (P  -  P  )_1V.  (73) 

From  Eq.  (30),  V  in  Table  3,  Pg  in  Table  4,  and  Eqs.  (73)  and  (68),  the  values 
of  Q  and  qg  in  Table  5  are  calculated.  Values  of  Vg  from  Eq.  (69)  are  also 
given.  By  setting  Pg  =  0  in  Eq.  (73),  the  charges  q.^  are  obtained  that  the 
conductors  would  have  at  the  same  potential  but  with  qg  =  0.  The  results  are  in 
Table  5.  Except  for  1^  =  1  pA,  the  magnitudes  of  q^  and  q ^  are  seen  to  be 
changed  no  more  than  9%.  So,  the  sheath  has  only  a  minimal  effect  on  the  con¬ 
ductor  charges.  By  using  Eq.  (6),  if  the  plasma,  2  and  3  were  absent  and  the 
potential  of  1  were  unchanged,  then  q^a>  shown  in  Table  5,  would  be  obtained 
for  the  charge.  The  results  show  that  the  presence  of  the  rear  section  causes 
about  89%  more  charge  to  accumulate  on  the  forward  section. 


4.8  Probe-Position  Potential 

If  q3  =  0,  then  V3  is  the  potential  of  the  position  of  the  probe  with  the 
probe  removed,  V  .  If  q_  *  0,  V  may  be  calculated.  A  simple  and  fairly 
accurate  method  is  to  subtract  from  the  part  of  the  potential  due  to  q^: 

Vppl  *  V3  -  P33V  <7,) 

Values  of  this  are  shown  in  Table  5.  The  proper  way  is  to  argue  that,  if  3  were 
not  there,  the  charges  would  be 

Q'  -  (P2  -  ps2)‘lv*  (75) 

where  q.',  V.  have  i.j  =  1,2,  and  P,  and  P  0  are  P  and  P  with  the  third  row 
and  column  deleted.  Then 

Vpp  =  (Pl3  *  P3s)ql’  +  (p23  '  P3s)q2’*  (76) 

Values  of  this  are  also  shown  in  Table  5.  They  are  seen  to  be  close  to  V  . , 

but  are  seen  to  differ  greatly  from  V^.  This  shows  that  |  q^  |  was  large  enough 

to  affect  Vj  greatly.  Except  for  Ifa  =  1  uA,  there  exists  a  value  of  conductance 

between  Sections  1  and  3  (less  for  the  ion  beams,  more  for  the  electron  beam 

than  !/z31)  for  which  q3  =  0  so  that  V3  =  Vpp.  This  may  be  seen  by  noting 

that  V  is  between  the  actual  V„  (Z_.)  and  V,  (Z„.  =«.)--  kT/e  for  ions 
pp  3  31  3  31 

except  for  I.  =  1  ^A)  and  V„  (Z.,,  =  0)  -  V  .‘or  electrons.  The  value  of  V 
r  b  3  Ji  I  pp 

for  the  I.  =  1  uA  case  should  be  0  or  slightly  negative;  -  5  V  is  unphysical. 

D  3 

Perhaps  this  indicates  that  1^  >  1  uA.  If  Katz  and  r.landell,  Section  4.1,  meant 

that  Vpp  is  40%  of  V^,  then  this  disagrees  with  the  above,  since  Vpp  =  10  to  15% 

of  Vj,  except  for  the  1^  =  1  yA  and  -10  ma  cases. 

4.9  Probe  Current 

In  the  absence  of  a  theory  for  a  probe  in  the  sheath  of  another  conductor, 
some  limiting  values  of  the  streaming  ion  or  electron  current  to  3,  I3>  are 
calculated.  Assuming  no  other  conductance,  this  current  is  balanced  by  the 
voltmeter  current 


z31  =  3  x  io10  a. 


Values  of  I. 


3vm  are  s*lown  *n  Table  5.  If  3  were  at  the  poten¬ 
tial  Vpp,  it  would  collect  approximately  a  fraction  of  1^  equal  to  the  fraction  of 
the  equipotential  surface  area  through  it  that  it  occupies.  Assuming  this  equipo- 
tential  surface  to  be  a  cylinder  concentric  with  the  forward  section,  this  area  is 


A  =  2nd  +  2nd2B 
c 


70.9  m 


(78) 


B  is  given  by  Eq.  (62),  where  a  =  a^,  b  =  b^  Rgc  =  d,  which  are  given  in 
Table  1.  If,  instead,  the  surface  were  assumed  to  be  a  sphere  concentric  with  1, 
the  area  would  be 


Ag  =  4n[d2  +  (b1  -  s) 2 ]  =  49.6  m2. 


(79) 


The  mean,  A  =  60  m  ,  will  be  used.  The  approximate  current  collected  if  = 


V  then  should  be 
PP 
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"a3  ^  =  1.9  x  10  Ix 


(80) 


Use  of  this  equation  yields  the  values  in  Table  5.  I3vm  is  seen  to  be  appreciably 
(generally  ~  60  times)  higher  than  i3Q  in  qualitative  agreement  with  being 
appreciably  less  (more  negative)  than  V  .  An  upper  bound,  is  obtained 

P  P  O  l » I 

the  same  way  but  substituting  the  maximum  impact  parameter 


p  =  a. 


eV, 


kT 


(81) 


in  place  of  a^.  Values  of  1-^,  also  shown  in  Table  5,  are  seen  to  have  far  too 

4 

large  a  magnitude.  Except  for  the  1^  =  -10  ^A  case,  this  is  mainly  because 
Eq.  (81)  gives  values  many  times  the  Debye  length,  whereas  only  particles  with 


an  impact  parameter  less  than  a  few  Debye  lengths  reach  t lie  probe.  Also,  this 
assumes  that  all  field  lines  are  radially  inward  to  3.  Since  this  is  incorrect, 
many  ions  with  impact  parameter  less  than  the  value  from  Eq.  (81)  will  miss  the 
probe.  The  values  of  1^  and  1^  for  1^  -  1  have  no  meaning  unless  the 
probe  is  inside  the  forward  section  sheath.  Assuming  the  probe  to  be  outside  of 


the  sheath  for  this  case,  spherical  Lam  theory  and  Eqs.  ( 38 ) - ( 46 )  give  i.  ->  X  +  1 


for  \'3  =  -30  V.  This  is  incorrect;  Lam  theory  must  be  invalid  since  r  <  1.  A 


fairly  accurate  value,  however,  may  be  obtained  from  a  careful  extrapolation  of 
Laframboise'st,  results  for  a  spherical  probe,  as  follows: 


25 


m 

as 


m 


Z', 


m 


■*v 


r  A 
c 


£1 


r.V? 


m 


53 


^  i  *  i 


v  v 


j&f 


Si 


yV 

-‘.V, 


:«aa 


(82) 


kx  =  1  +  .481C  -  .108C2 

k2  =  .0331C  +  ,078C2 

k  -  1 

-X  =  kx(i  -  1)  -  —  ln[  1  +  k2(i  -  1)]. 


(83) 

(47) 


The  four  numbers  in  Eqs.  (82)  and  (83)  are  chosen  so  that  Eq.  (47)  agrees  with 
Laframboise  for  C  =  -2  and  .3  at  x  =  25,  and  for  C  =  .5  at  x  =  20  and  25.  As¬ 
suming  that  Laframboise's  numbers  are  correct,  for  0  <  c  <  .5  and  0  <  x  £  25, 
these  equations  yield  x  correctly  probably  within  ,2%.  For  0  <  x  s  25,  the  error 

is  certainly  less  than  .8%  for  C  =  .2,  1.9%  for  C  =  .3,  and  .4%  for  C  =  .5.  As- 
9  -3 

suming  n  =10  m  ,  C  =  .372  and  these  equations  give  I~  =  16  nA.  Assuming 

°  8  -3  * 

n  =  5  x  10  m  (at  the  sheath  edge),  they  give  I,  =  8.4  nA.  This  is  still  an 

order  of  magnitude  larger  than  I3  .  The  contribution  to  1^  of  secondary  and 
scattered  electrons  from  the  forward  section  is  negligible.  Perhaps  this  dis¬ 
crepancy  is  due  to  neglecting  the  effect  of  the  magnetic  field  that  decreases  the 
return  current.  Alternatively,  perhaps  it  indicates  that  the  probe  was  inside  the 
sheath  where  these  probe  theories  are  invalid.  This  would  favor  the  larger  value 
of  I.,  for  which  the  probe  would  more  likely  be  inside  the  sheath.  Although 
unlikely,  it  is  possible  that  the  extrapolation  yields  much  too  high  a  value  of 
current . 

As  |lbi  is  increased  up  to  the  minimum  value  for  saturation,  jVg|  increases. 
As  |Ib|  is  increased  further,  | Vg |  decreases  and  V3  asymptotically  approaches 
~  SkT/e  as  |lb|  -*■  «,  where  S  =  1  for  an  ion  beam,  and  -1  for  an  electron  beam. 
The  reason  for  |Vj|  decreasing  is  that,  above  minimum  saturation  as  |lb|  in¬ 
creases,  if  V3  remained  constant,  Igvm  and  the  fraction  of  return  current  col¬ 
lected  by  the  probe  would  remain  constant,  so  the  amount  collected  would  inci’ease 
and  the  two  currents  would  not  balance.  Thus,  Vg  must  change  to  decrease  the 
fraction  collected,  and  to  increase  I„  to  maintain  balance.  So,  for  a  2  keV  ion 
beam,  |V3|  would  be  maximum  at  Ib  *  26  yA;  for  a  90  eV  electron  beam,  V3 
would  be  maximum  at  =  -853  yA.  Since  the  electron-beam  case  has  Vr3  =  0,  it 
is  seen  to  be  much  more  highly  saturated  than  the  Ib  =  374  yA  case. 

4.10  Charging  Time  Constant 

Now  the  charging  time  due  to  current  and  capacitance  is  calculated.  Lower¬ 
case  letters  are  used  for  instantaneous  values  of  variables,  and  upper-case  letters 
for  steady-state  values.  The  effect  of  the  probe  is  neglected.  Differentiating 


26 


iulwij  &JHVIU  imu 


■WOTKmvansiCHjrvijnvianancnvunwvwvnn 


the  first  two  rows  of  Eq.  (2)  with  respect  to  time,  and  equating  to  the  current: 
^1  =  ClCl  +  c12°2  -  (cll  +  c12);'ls  =  il~  lb  + 


v2  ‘  V1 


'21 


(84) 


q2  =  c^Vj  +  c2‘,^2  ~  (C12  +  C22)'r 


=  - 


V2  V1 


Is  ‘2 


(85) 


The  assumptions  are  made  that: 


_  +  <*2 


'Is 


Q~  *  Is 
s 


Q1  +  q2 


Is 


(86) 


V1  -  V10  r,  V2  -  V1  1 

=  [ b "  J 


(87) 


= 


V2  ~  V20 


V  -  V 
2  V1 


v  ^~V 
v2  v20 


'21 


(88) 


Approximate  values  of  the  conductor  potentials  before  the  beam  is  turned  on, 
V.q,  need  to  be  known.  For  them,  Eqs.  (38)  to  (46)  are  used  with  I  =  0  (or 


(V20  -  W' 


Interpolating  numbers  in  Laframboise  Table  6c,  V.n  =  -.214 


=  -.196  V  is  obtained. 


10 

The  beam 


V  is  obtained.  Using  Eqs.  (47)  to  (49),  V£q 
is  assumed  to  be  turned  on  at  t  =  0  and  to  be  a  step  function;  that  is,  to  have 
a  rise  time  much  shorter  than  the  charging  times.  The  solution  is  found  to  be 


q.  =  Q.  +  A.e  t/xl  +  B.e  t/x2. 
n  l  l 


(89) 


The  algebra  is  lengthy,  and  the  expression  for  (functions  of  the  constants  in 
Eqs.  (84)  to  (88),  the  solution  of  a  quadratic  equation)  is  long.  One  time  con¬ 


stant  is  10%  (1  or  2%  for  1.  =  1  pA  and  20%  for  I.  =  -10  uA)  of  the  other. 

b  b  3-3 

Values  of  the  larger  time  constant  are  shown  in  Table  6.  For  nQ  =  10  cm 


the 


ion  and  electron  plasma  periods  are  .76  ms  and  3.5  us,  respectively.  The 
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Table  6.  Charging  Time  Constant  and  Initial  Ion  Transit  Time  Through  the  Sheath 


11  -  — 1 

*b 

1 

m 

1 

8 

a 

.O 

00 

8c 

9 

10.9 

12 

374 

1 

I— * 

o 

uA 

T 

8.0 

2.6 

6.0 

5.4 

5.7 

5.4 

5.6 

5.9 

8.8 

.040 

ms 

tts(T) 

.90 

2.2 

2.4 

(9.4) 

(.5) 

ms 

w 

(13.8) 

(1.4) 

ms 

ttc(T) 

.93 

1.7 

2.1 

2.4 

.19 

ms 

w 

.32 

ms 

charging  time  constant  is  seen  to  be  much  larger.  A  simple,  approximate  way  of 

4 

calculating  t  is  given  in  Reference  2.  Except  for  I.  =  374  and  -10  uA,  it  gives 
values  within  10%  of  those  in  Table  6.  For  Ib  =  374  and  -10  pA,  the  values  in 
Reference  2  are  not  values  of  t,  which  should  have  been  7.6  and  .0078  ms, 
respectively,  but  are  close  to  the  values,  .65  (assuming  6  =  .18)  and  .0012  ms, 
respectively,  (from  the  simple  way)  of  tj,  the  time  for  the  potential  of  the 
forward  section  first  to  reach  Vj  =  Eb.  (Vj  then  overshoots  Vj  before  finally 
settling  at  V^)  Values  of  t^  also  could  be  calculated  easily  by  the  above  more 
accurate  method.  At  least  with  the  simple  way,  x  is  inversely  proportional  to 
(1  +  6)Ij,  so  it  has  this  dependence  for  such  saturated  cases;  however,  since 
(1  +  6)1^  =  1^  -  Ij  for  unsaturated  cases,  t  is  then  independent  of  6,  given  Ib> 
Actually,  dl^/df-V^  for  ion  beams  decreases  with  increasing  (-Vj);  similarly, 
dij/dt-Vj)  probably  decreases  with  increasing  (-v).  Taking  this  into  account 
would  lead  to  the  charging  being  nonexpotential.  The  treatment  here,  however, 
should  give  approximately  the  correct  behavior. 


4.11  Slieatli  Ion  Initial  Transit  Time: 


The  approximate  traverse  time,  tt>  is  calculated  for  the  first  ions  after  the 
beam  is  turned  on  to  travel  radially  between  Rg  and  R^,  the  assumed  forward  sec¬ 
tion  radius.  In  Reference  2,  the  sheath  was  assumed  spherical,  and  of  radius 
equal  to  the  value  tabulated.  It  was  concentric  to  a  spherical  conductor  of 

.5  m  radius  representing  the  forward  section.  The  values  obtained  of  t  =  t 

^  i  is 

for  t'  =  x  or  t.  are  shown  in  Table  6.  Those  for  I,  =  374  and  -10  uA  should 
ts  b 

be  discounted.  They  are  much  too  large  because  the  values  of  Rg  are  too  large. 

Now  the  sheath  is  assumed  cylindrical  and  concentric  with  the  actual  forward 

section.  Figure  2.  For  the  ion-beam  cases,  the  ion  is  at  R  =  Rsc  with  R  =  -v^ 

at  t  =  0.  Eq.  (A28),  Appendix  C,  is  integrated  numerically  to  t  =  tf ( t ' )  where  R 

=  a^.  For  the  electron-beam  case,  the  ion  is  at  R  =  a^  and  R  =  0  at  t  =  0;  Eq. 

(A28)  is  integrated  numerically  to  t^  ( t  1 )  where  R  =  Rgc.  Some  results  for  t  = 

t..  and  t'  =  t  or  t.  are  shown  in  Table  6.  The  values  of  I,  =  8,  and  9  yA  are 
tc  tc  b  b 

expected  to  be  nearly  identical.  The  results  for  the  two  models  are  seen  to  be 

about  the  same.  The  values  for  t^  for  a  spherical  sheath  around  the  cylindrical 

forward  section.  Figure  3,  would  be  a  little  larger  since  Rgs  -  a^  is  larger  than 

R  -a,,  t.  is  seen  to  be  smaller  than  t  except  for  the  I.  =  -10  ma  case.  So, 
sc  1  t  ^  b 

t  is  the  minimum  time  in  the  ion-beam  cases,  and  .32  ms  in  the  electron-beam 
case,  for  steady  state  to  be  reached. 

5.  CONCLUSIONS 

The  charging  data  are  analyzed  for  eight  positive  ion  cases,  and  one  elec¬ 
tron-beam  emission  case,  from  a  three-conductor  rocket  in  the  lower  F-region  at 
night.  Calculation  of  sheath  radii  show  that  the  spherical  emissive  probe  was  in 


the  sheath  of  the  forward  section,  in  agreement  with  Katz  and  Mandell's  interpre- 
3 

tation,  although  the  1^  =  1  yA  beam  current  case  is  a  possible,  but  unlikely, 
exception.  How  potential  coefficients  for  any  number  of  conductors  can  be  mea¬ 
sured  is  indicated.  Since  the  coefficients  were  not  measured  for  the  rocket,  and 
the  geometry  was  sufficiently  simple,  values  of  them  and  of  sheath  potential 
coefficients,  both  dependent  only  on  the  geometry,  are  calculated  fairly  accur¬ 
ately.  An  empirical  formula  found  for  the  potential  due  to  a  charged  conducting 
cylinder,  Smythe's  charge  distribution  on  a  conducting  cylinder,  his  empirical 
formula  for  the  capacitance,  and  three  gedanken  experiments  are  used  to  deter¬ 
mine  these  coefficients. 
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By  using  these  coefficients,  understanding  was  gained;  this  included  the  charge 
on  the  conductors  and  in  the  sheath,  the  potential  at  the  position  of  the  probe 


without  the  probe,  V  ,  and  the  charging  time  constant. 


An  empirical  formula  accurate  to  .2  to  4%  in  different  ranges  of  probe  radius 
to  Debye  length  ratio,  i ,  up  to  5  is  found  for  Laframboise's  spherical  probe  re¬ 


sults.  Empirical  formulae  accurate  to  1"  over  the  whole  range  are  found  for  F 


,2  . 


and  G.  F  is  useful  in  Lam's  theory,  and  in  lieu  of  (-  a)  in  Langmuir-Blodgett’s 

,2  . 


spherical  probe  theory.  G  is  useful  in  Lam’s  theory,  and  in  lieu  of  (-6)  in 
Langmuir-Blodgett's  cylindrical  probe  theory. 

The  spherical  probe  theories  of  Laframboise  and  Lam,  and  the  rear  sec¬ 
tion  voltmeter  current,  lead  to  the  potential,  V of  the  rear  section.  The  values 
are  negative  and,  since  the  contributions  from  the  various  charges  nearly  cancel, 
are  of  small  magnitude  in  qualitative  agreement  with  a  relevant  calculated  value  by 
Katz  and  Man  dell.  In  general,  however,  the  magnitudes  are  appreciably  larger  due 
to  voltmeter  conductance.  These  and  measurements  of  give  absolute  values  of 
potential  of  the  forward  section,  V^,  and  probe,  V^.  For  ion-beam  cases,  the  front 
part  of  the  rear  section  is  concluded  to  have  a  positive  surface  charge  density, 
and  the  rear  part  a  negative  surface  charge  density.  The  latter  disagrees  with 


Katz  and  Mandell.  Values  of  return  current  are  calculated  from  V1  by  using  Lam's 

1  4 

probe  theory  for  a  cylinder.  Except  for  the  two  high-beam  current  (374  and  -10 


ijA)  cases,  they  agree  roughly  with  Ib  -  1^.  So,  probe  theory  accounts  for  the 
high  value  of  |Vj|  for  these  unsaturated  cases  (e  ]  V1  ]  <  Eb,  the  particle  energy 
of  the  beam).  It  also  shows  that  the  two  high-beam  current  cases,  especially  the 
electron  case,  are  highly  saturated  (Ib  much  greater  than  necessary  for  ejVjJ  = 


Eb>;  most  of  the  beam  current  is  reflected  near  the  sheath  edge  back  to  the  for¬ 


ward  section. 

For  the  ion-beam  cases,  although  all  secondary  electrons  from  the  forward 
section  escape,  it  appears  that  they  have  a  negligible  effect,  other  than  de¬ 
creasing  the  return  current  and  sheath  radius  slightly.  For  the  electron-beam 
case,  many  secondary  and  scattered  electrons  come  from  the  forward  section  sur¬ 
face,  but  they  cannot  escape,  and  have  very  little  effect,  even  on  the  sheath 
charge  density.  The  effect  of  ionization  by  ions  is  negligible,  that  by  electrons 
is  expected  to  be  small  to  moderate,  and  the  effect  of  magnetic  field  on  electrons 
is  expected  to  be  quite  appreciable. 


30 


m 


XI 


Si 


t  A 


vN 


£ 


A.’M 


Values  of  charge  on  the  rocket  conductors  and  in  the  sheath  are  calculated. 
They  are  not  greatly  changed  by  the  sheath  charge.  The  charge  on  the  rear  sec¬ 
tion  must  always  have  the  opposite  sign  and  a  smaller  magnitude  than  that  on  the 
forward  section.  The  presence  of  the  rear  section  is  found  to  increase  the 
charge  on  the  forward  section  by  about  89?.  The  potential  of  a  conductor  is  con¬ 
cluded  to  be  determined  mainly  by  the  current  from  it,  and  the  charge  on  it  to  be 
determined  by  potentials  and  induction  (coupling).  Both  a  simple  and  an  accurate 
way  are  used  to  calculate  V  .  It  is  much  higher  (less  negative)  than  V^.  This 
is  due  to  cjg,  the  charge  on  the  probe  (that  is  related  to  the  equilibrium  probe 
current).  V  is  10?,  to  15?  of  Vp  rather  than  40?  as  apparently  interpreted  by 
Katz  and  Mandell.  For  each  ion-beam  case,  except  1  yA,  a  value  of  probe  voltmeter 
resistance  exists  for  which  would  equal  0,  and  V  would  equal  V^.  Upper  and 
lower  limits  of  probe  current  from  the  plasma  are  calculated  that  bracket  the  bal¬ 
ancing  probe  voltmeter  current,  1^  ,  For  1^  =  1  yA,  the  value  of  return  current 

to  the  probe  calculated  by  probe  theory  is  an  order  of  magnitude  larger  than  I,  . 
This  is  probably  due  to  neglecting  the  effect  of  the  magnetic  field,  but  may  also 
be  due  to  plasma  current  to  the  forward  section,  and  electric  field  from  the  for¬ 
ward  section,  i.e.,  to  the  probe  being  inside  of  the  forward  section  sheath.  As 
the  beam  current  increases  beyond  minimum  saturation,  approaches  zero  and, 
then,  at  extreme  values,  kT /e  beyond  zero.  The  electron-beam  case  is  in  accord¬ 
ance  with  this;  further  verification  would  be  desirable. 

The  charging  time  constant  due  to  capacitance  is  an  order  of  magnitude  larger 
than  the  plasma  period,  and,  except  for  the  electron-beam  case,  is  larger  than  the 
initi;il  ion  transit  time  through  the  steady-state  sheath.  So,  the  minimum  time  to 
reach  steady  state  is  determined  by  the  charging  time  due  to  capacitance  for  an 


ion  beam,  and  by  the  ion  transit  time,  a  third  of  a  ms,  for  the  electron-beam 


case.  The  time  for  the  374  or  -10  yA  case  to  become  saturated  was  loss  than  a 


plasma  period. 

In  spite  of  many  approximations,  this  analysis  is  believed  to  be  qualitatively 
and  approximately  quantitatively  correct,  at  least  for  the  ion-beam  cases.  Since 


"  <.  10  for  the  cases  here,  a  better  probe  theory  than  Lam's,  perhaps  that  of 


Bernstein  and  Rabinowitz.  '  should  In'  used  for  potentials  out  of  range  of  those 


of  Laframboise. 


13.  Bernstein,  I.  B.,  and  Rabinowitz,  1.  (1959)  Theory  of  electrostatic  probes 
in  a  low  density  plasma,  I’hys.  Fluids  2:112-121. 


For  the  electron-beam  case,  a  treatment  is  needed  that  includes  the  effect 

of  magnetic  field,  and  probably  ionization,  for  reliable  results.  An  accurate 

(numerical)  solution  of  Poisson's  equation  would  make  an  interesting  and  perhaps 

valuable  comparison  with  the  results  in  this  report.  The  analysis  techniques  used 

here,  applied  to  a  greater  variety  of  data,  should  produce  much  more  understanding 

12 

From  the  results  in  this  report,  and  those  of  Lam-Jaycor,  one  may  induce 
the  following  general  steady-state  properties  of  vehicles  in  the  ionosphere  at 
potentials  small  enough  not  to  cause  breakdown  or  appreciable  ionization  of  back¬ 
ground  neutrals.  For  a  given  beam  current  and  shape  of  vehicle  with  an  all¬ 
conducting  surface,  as  the  size  increases,  the  sheath  thickness,  vehicle  potential, 
and  charge  decrease.  For  sheaths  that  are  thick  compared  to  vehicle  size,  the 
charge  decreases  very  little.  For  a  given  beam  current  and  vehicle,  as  the  frac¬ 
tion  of  the  surface  that  is  covered  with  dielectric  is  increased,  the  potential  and 
charge  on  the  conducting  part  increase. 

It  is  hoped  that  this  report  will  provide  useful  considerations  and  techniques, 
and/or  stimulate  ideas  leading  to  better  ones  for  experiments  and  analyses  of 
object-charging  in  laboratory  and  space  plasmas. 
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Appendix  A 

Potential  Due  to  a  Charged  Conducting  Cylinder 


The  following  is  an  empirical  formula  for  the  potential  of  a  point  R,  z  (see 
Figure  2)  due  to  a  conducting  right,  circular,  cylinder  of  radius  a  with  a  charge 
q,  b  >  a,  R  >  a: 


V  =  g^?  fl 


R*  1  + 


2b  +  R  z 


a  in  meter,  q  in  Coulomb,  V  in  Volt. 


When  b  <<  |z|  and  b^  -  Rjzl, 


V  =  g 


R  \  1  + 


2b  +  R  z 


When  (b  <<  Izl  and  b  <  R)  or  b  R, 


-  1  . 


When  -g  «  In  ^  «  —  and  |z|  <  b. 


When  z  =  0,  R  =  a,  k  =  1.777,  and  a/2  <  b  <  8a,  then  q/V  is  within  1.25%  of  the 
correct  capacitance. 

Pl03  1  *s  obtained  from  Eq.  (Al)  with  a  =  .19  m,  k  =  1.777,  b  =  1.27  m 
R  =  1.71  m,  z  =  1.01  m  giving 


p103  ~  •2605  nf‘ 


p123  1  *s  obtained  from  Eq.  (Al)  with  a  =  .19  m,  k  =  1.777,  b  =  1.475  m, 
R  =  1.71  m,  |  z  |  =  .805  m,  giving 


=  .261-  nf. 


(A6) 
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Appendix  B 


Potential  Due  to  the  Sheath 


B.l.  Cylindrical  Sheath 

The  potential  at  R  on  the  median  plane  z  =  0  due  to  a  cylindrical  sheath 
(Figure  2)  is  calculated.  From  Eq.  (Al),  the  potential  due  to  an  infinitesimal 
cylindrical  shell  at  R'  <  R  and  z'  is 


dV  = 
< 


_kdqs  [7.  +  Z' 

g“HT’L\  * ) 


where 


!’  =  R'[1*  (!  - ')  (fr)  b'a]  ■ 


In  this  section,  a  =  a^  and  b  =  b^. 

The  potential  due  to  a  shell  at  It'  >  R  is  assumed  to  be 


This  implies  that  the  charge  distribution  is  that  for  a  conducting  shell.  The 
sheath  volume  is 


Vol  =  wR22B  -  i a2 2b, 
8 


(A10) 


where  B  is  given  by  Eq.  (62). 


=  =  2"p[3  +  (2J  ■3)(^r)6^i]R'2dR,• 


(All) 


where 


2trp[R2B  -  a2b]  =  q  . 
s  s 


( A12) 


This  assumes  that  the  average  charge  density,  p,  is  the  same  in  each  shell 

volume.  Errors  due  to  this  assumption  are  believed  to  be  small  since  Lafram- 
6 

boise,  Figure  33,  shows  p  to  be  nearly  constant  in  the  sheath  for  £  =  2 
(±  13%  to  4Xp)  and  5  (±9%  to  8x^  from  the  probe  surface)  and  since  and 
are  not  greatly  changed  by  neglecting  the  plasma  charge  entirely.  The  potential 
due  to  each  shell  is  of  the  form 


dV<  =  2*gkpu(R\R)dR' 


(A13) 


dV>  =  2wgkpw(R')dR'. 


( A14) 


Empirical  formulae  are  found  for  u  and  w : 


u(R',1.71)  =  .008  +  .47R’  +  .648R’ 


( A15) 


w(R')  =  .247  +  1.427R'. 


(A16) 


Eq.  (A15)  is  correct  within  1%  for  .19  <  R'  <  1.71  m.  Eq.  (A16)  is  correct 
within  9%  for  .19  <  R'  <,  .3  m  and  within  2%  for  .3  m  <  R\ 

The  potential  per  unit  sheath  charge  at  the  side  of  the  rocket  due  only  to 
sheath  charge  is  taken  to  be 


^ -•  -s  v  -v  - 


( A17) 


p,  =  /  dV  /q  . 
Flsc  >  Ms 


The  potential  per  unit  sheath  charge  at  the  probe  due  only  to  sheath  charge  is 
taken  to  be 


P3sc  =/  dV</qs  +  f  dV>/qs 

JSC  a  S  1.71 


( A18) 


for  the  probe  inside  of  the  sheath  and 


p3sc  =  f  dV</qs 
a 


( A19) 


for  the  probe  outside  of  the  sheath  (only  used  for  1^  =  lm  yA). 


B.2.  Spherical  Sheath 


As  R  increases,  the  sheath  boundary  becomes  less  cylindrical  and  more 
s 


spherical.  For  R  >>  b,  this  boundary  (Figure  3)  is  nearly  spherical,  so  is 
s 


assumed  spherical.  The  potential  per  unit  sheath  charge  at  the  side  of  the 
rocket  due  only  to  sheath  charge  is  taken  to  be 


plss  Vss^a^qs, 


(A20) 


and  that  at  the  probe  is  taken  to  be 


=  V,s<1-7»'v 


(A21) 


where  V  (R)  is  derived  on  pp.  34  and  35  of  Reference  2  and  given  in  Eqs. 
ss 


( A35)-(A37)  there. 


iV  vM*.  ^  A  <,  <■.  /.v. 


m. 


& 
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Appendix  C 

Initial  Ion  Transit  Time  Through  the  Sheath 


First,  the  potential,  V(R),  is  calculated  in  the  sheath  at  a  radius  R,  along  a 
radial  line  passing  through  the  center  of  the  forward  section  (z  =  0).  It  is  the 
sum  of  the  potential  at  R  due  to  qg,  q^,  and  q2  (the  effect  of  is  small  and  is 
ignored) : 


R  R 

V(R)  =  27.gkp[y*  u(R',R)dR'  +  f  SCw(R')dR'J 
a  R 


+ 


g  1 

a 


R 

a 

17  bAroi  1 

kq2  \( 

lH  iJ 

*  [{' 

f'-l 


(A22) 


from  Eqs.  (Al),  (A13),  and  (A14),  where 


S=  1  + 


<bl  +  b2> 

2b*  +  R2(bj  +  b2)2 


(  A23) 


The  electric  field  is  assumed  to  increase  everywhere  exponentially  with  time,  t, 
after  beam  turn  on  with  a  time  constant  x\  So 
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where 

R3u 

KR)  =f  1^-dR’, 

fl 

and 

u(R,R)  -  w(R)  =  0. 

The  empirical  formula 

I(R)  =  -A(R  -  a)  -  B  In  R  +  C(I  -  , 

A  =  .4347,  B  =  .02784,  C  =  .02093 

is  accurate  to  .3%.  The  equation  of  motion  is 
R  E(R,t), 


( A24) 


( A25) 


( A26) 


(A27) 


(  A28) 


where  e  and  M  are  the  magnitude  of  charge  and  the  mass  respectively  of  the  ion, 
assumed  to  be  NO+. 


